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Abstract. OB-stars have the highest luminosities and strongest stellar winds of all stars, which enables them to interact 
strongly with their surrounding ISM, thus creating bow shocks. These offer us an ideal opportunity to learn more about the 
ISM. They were first detected and analysed around runaway OB-stars using the IRAS allsky survey by van Buren et al. ( 1995 I. 
Using the geometry of such bow shocks information concerning the ISM density and its fluctuations can be gained from such 
infrared observations. As to help to improve the bow shock models, additional observations at other wavelengths, e.g. Ha, are 
most welcome. However due to their low velocity these bow shocks have a size of ~ 1°, and could only be observed as a whole 
with great difficulties. In the light of the new Ho- allsky surveys (SHASSA/VTSS) this is no problem any more. 
We developed different methods to detect bow shocks, e.g. the improved determination of their symmetry axis with radial 
distance profiles. Using two Ha-allsky surveys (SHASSA/VTSS), we searched for bow shocks and compared the different 
methods. From our sample we conclude, that the correlation between the direction of both proper motion and the symmetry 
axis determined with radial distance profile is the most promising detection method. 

We found eight bow shocks around HD 17505, HD 244 30, HP 48099, HD 57061, HD 92206, HD 135240, HD 149757, and HD 
158186 from 37 candidates taken from van Buren et al. ( 1995 1. Additionally to the traditional determination of ISM parameters 
using the standoff distance of the bow shock, another approach was chosen, using the thickness of the bow-shock layer. Both 
methods lead to the same results, yielding densities (~ 1 cnT 3 ) and the maximal temperatures (~ 10 4 K), that fit well to the 
up-to-date picture of the Warm Ionised Medium. 
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1. Introduction 

OB-stars are the most massive and luminous stars known with 
masses greater than 10M o and effective temperatures rang- 
ing from 10000K up to 50 000 K. During their short lifetime 
(<20xl0 6 yr) OB-stars lose mass at rates of M » 10~ 7 - 
1O~ 5 M yr _1 (Lamers & Cassinelli [1989 1. This stellar wind 
with velocities of » 1 000 - 3 000 km s _1 transfers a great 
amount of mechanical energy to the surrounding ISM, com- 
parable to a supernova explosion. As a result, OB-stars cre- 
ate a stellar bubble (Castor et al. I1975> which structures the 
Interstellar Medium (ISM). This spherical bubble is altered 
when the OB-star is in motion, as wind and ISM interact di- 
rectly. The resulting nebula is known as a bow-shock nebula. 
Recently many new bow shocks (size ~1') created by fast mov- 
ing neutron stars or pulsars have been detected and analysed. 
Their geometry could be successfully used to gain informa- 
tion about the ISM, such as density and temperature (Gaensler 



et al. 2001 1. Thus the theoretical models of such bow shocks 
have become quite exact recently, enabling us to glance at the 
missing link between density fluctuations seen in H i (at scales 
~ 0.1 - 200pc) and towards pulsars (at scales ~ 5 — 100 AU). 
Bow-shock nebula around OB-stars can also be used as ISM 
probes, if observable or rather detectable, as described in the 
following. 

In 1988 van Buren & McCray detected structures around OB- 
stars and Wolf-Rayet stars using the 60 /mi allsky survey of 
the InfraRed Astronomical Satellite (IRAS). These images re- 
vealed an arc-like structure and a high colour temperature, pos- 
sibly being bow shocks. A more complete sample of 188 run- 
away OB-stars was analysed by van Buren et al. ( 1995) (here- 
after VB) using the IRAS allsky survey which lead to the detec- 
tion of 58 bow shocks. Because of these numerous detections 
and the alignment of the symmetry axis of the structures along 
the direction of proper motion of the central star, they could 
only be bow-shock nebula. 
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Fig. 1. Structure of a stellar bubble created by a stationary star 
(left) and a star in motion (right). The different regions A, B, 
C, and D and the geometry are described in Sect. 13 

Due to the low resolution of IRAS ( — 1 '5) it is difficult to de- 
termine the exact location of the bow shock and its symmetry 
axis. Furthermore, VB could only use the Hipparcos input cat- 
alogue (HIC) to determine the proper-motion direction of the 
central star. As the Hipparcos catalogue is now completed, its 
astrometrical data together with the almost complete Ha all- 
sky surveys Virginia Tech Spectral Survey 1 (VTSS; Dennison 
et al. 119971 and Southern Hemispheric Ha Sky Survey Atlas 1 
(SHASSA; Gaustad et al. l2001i . the VB sample is reanalysed 
here, now using the Ha emission line. 

The development from a stellar bubble to a bow shock is il- 
lustrated in Sect. 2. The used Ha-data and its acquisition is 
described in Sect. 3. Sect. 4 explains the analysis of the bow 
shocks and the improved methods developed to detect them. In 
Sect. 5 the results of the observations are given and in Sect. 6 
ISM parameters derived, followed by a discussion in Sect. 7. 
Finally, Sect. 8 draws a conclusion concerning the main points 
of this paper. 

2. Scenario 

For a better understanding of the methods used, the scenario of 
a bow shock and the source of the high velocity of the OB-stars 
are given. 

The evolving stellar bubble, as illustrated in Fig. [2 is in its 
longest snowplow stage divided into four parts centred upon 
the OB-star (Castor et al. I1975> : The innermost area (A) is 
an unshocked and freely expanding stellar wind, followed by a 
larger area (B) of shocked stellar wind. These hot regions, as 
if they were a snowplow (hence the name of this stage), have 
pushed together a thinner area (C) of shocked ISM. All regions 
are embedded in the unshocked ISM (D). 
The models describing such stellar bubbles assume that the 
OB-star is stationary with respect to the ISM. However as all 
stars have a proper motion, so do OB-stars. A special pop- 
ulation of OB-stars is known as runaway OB-stars (Blaauw 
I196U . They are defined as having a proper motion greater than 
30kms~'. This criteria was chosen to discern them from non- 
runaway OB-stars, which have a velocity dispersion of about 
lOkms- 1 . 

As runaway stars are often found in isolated regions, their high 
1 Supported by the National Science Foundation 



velocity cannot be explained as a motion within a stellar clus- 
ter. Two scenarios explaining observed properties of runaway 
OB-stars are favoured: 

Firstly, the Binary-Supernova-Scenario (BSS) as described 
by Blaauw (1961 1. The partner of the OB-star explodes as 
a Supernova (SN). Thereby the OB-star is set free with 
its typical orbital velocity of 30-150kms _1 . And secondly, 
the Dynamical-Ejection-Scenario (DES) proposed by Hoffer 
(1983). In this Scenario the collision of two binary systems 
leads to the ejection of one star with a velocity of up to 
200 km s" 1 . 

Taking the motion of the star (V„) into account, regions A, B, 
and C will still be spherical as long as the stellar velocity is 
smaller than the sound velocity within B. However B and C are 
no longer centred on the star, as shown in Fig ^ If region A 
and D do not interact directly, this would be the only alteration 
to the model of the stellar bubble. But as soon as the star enters 
denser regions of the ISM, like molecular clouds, the cooling 
of region B becomes more effective. This leads to a collapse of 
B and C within timescales smaller than the lifetime of the OB- 
star, and the approach of region A and D. As molecular clouds 
are not frequently encountered, region A and D can only inter- 
act when the offset of A to B and C is Vj « R B (t)-R A (t). Where 
the timescale is dependent upon the velocity of the OB-star and 
the density of the surrounding ISM. Taking typical values of 
the lifetime and velocity V,,=39kms~ 1 this leads to a density 
constraint of n > 0.017 cm 3 in which A and D can interact 
directly. 

In the case of directly interacting unshocked stellar wind 
and unshocked ISM, the geometry is changed completely. As 
Wilkin (1996 1 describes, the ram pressure of both media can 
be balanced directly and result in a bow shock. This bow shock 
is axi-symmetric along the direction of proper motion and can 
be approximated by a parabola. The two layers B and C of the 
model above are mixed due to turbulence and plasma instabil- 
ities leading to a single layer in which the material of the ISM 
and stellar wind moves along the bow shock. The material in 
this has experienced a nearly isothermal shock, so its density 
is higher than that of the surrounding ISM. This leads to the 
creation of warm interstellar dust best seen in 60 /mi, and the 
OB-star in the centre leads to the ionisation of the layer emit- 
ting Ha. 

3. Data 

3.1. Selection 

The data used for this program were taken from the SHASSA 
and incomplete VTSS allsky surveys. SHASSA contains the 
southern hemisphere up to £= 15° and VTSS the northern hemi- 
sphere down to 6 =-15°. Both surveys were made with a CCD 
detector and a fast photo-objective of ~55 mm at ~f/1.4, lead- 
ing to a field of view of ~13°. All images were integrated 
~25 min resulting in a detection limit down to ~0.75 rayleigh 
(1 R = M Ph°t™s \ £) ue t0 t h e different pixel sizes of the detec- 

V In cm" s sr / r 

tors SHASSA has a resolution of 0^8 and VTSS of 1 .'6. 

The Ha-sample used consists of the O-stars taken from the 

VB sample with data from either VTSS or SHASSA, ensur- 
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Fig. 2. Left: IRAS 60 fjm excess map of HD 1 35240. Right: Overlay of Ha image in grey-scale and the IRAS 60 jim excess map 
in white contours. Both are shown with their slanted equatorial-coordinate system (J2000) and inverted grey-scale. 



ing a sufficient Lyman continuum flux to ionise the bow-shock 
layer. 

We searched for bow shocks around these 37 candidates of the 
Ha-sample within the SHASSA and VTSS Ha survey. Due 
to always present background nebulosity, it had to be ensured 
that the structures seen were really bow shocks. As described 
at the end of Sect. [5] the bow-shock layer should be visible 
in 60 /mi as well as in the Ho- emission line. Therefore, the 
Ha images were compared with the 60 //m IRAS images of the 
same region. Though the IRAS images show a great amount 
of nebulous emission, the nebulosity can successfully be sub- 
tracted using the 100 /mi images of IRAS. This correction has 
been done based on the recipe of VB with the creation of IRAS 
60 /mi excess maps, shown for the example of HD 135240 in 
Fig.0 

To compare both images with each other, the contours of the 
IRAS 60 /mi excess map were overlaid upon the Ha images 
(see Fig.|2J. We used this image to decide whether a bow shock 
seen in the IRAS 60 fmx excess map is also present within the 
Ha image. The eight bow shock detections and a short descrip- 
tion of the comparison using the overlays are given in Table 
H The IRAS 60 /mi excess images and their corresponding Ha 
images are also shown in Fig.[S]and Fig- El 
To check for the positioning quality of the Ha image com- 
pared to the IRAS 60 //m excess maps in the overlays, we 
could not use the nebulosities as criteria. As we are using 
these overlays to see coinciding positions of nebulosities within 
both images this would be misleading. Better criteria are point 
sources like stars which have to be visible with IRAS and in 
VTSS/SHASSA, thus, we used the positions of M giants as a 
reference. Using this method, a deviation of the positional off- 
sets of cr = 20" was measured, which is sufficiently smaller 
than the resolution of both images. 
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Fig. 3. Median filtered section of the Ha image containing the 
bow shock around HD 135240. The coordinate system and 
grey-scale as in Fig. [2] 



To analyse the Ha images of the selected bow shocks only the 
interesting region was extracted. As only the search for bow 
shocks using the overlays requires the best possible resolution 
of 1.6' (VTSS) and 0.8' (SHASSA), the images could for fur- 
ther analysis be median filtered with a 5 x 5 pixel wide box (see 
Fig-Efor the example of HD 135240). The resulting improve- 
ment of S/N leads to a decreased resolution of 8' (VTSS) and 
4' (SHASSA). 

3.2. Distances 

To convert angular sizes into linear sizes, the distances 
of the eight bow-shock candidates had to be determined. 
Additionally, the interstellar absorption had to be calculated to 
fit the brightness profile correctly. 
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Table 2. Photometry and astrometry of the central stars of the bow-shock candidates, oig and my are the apparent B and V 
magnitudes, /u a and /j t ; give the proper motion along the rightascension and declination axis, and V r is the radial velocity. The 
spatial velocity V, the symmetry axis position Q a and its inclination La are derived from the astrometric data. The distance r is 
derived from the spectral parallax. 



star 


SC 


multiplicity 


m B 


my 


Ha 




V r 


V 


r 


e a 










[mag] 


[mag] 


[mas yr~'] 


[mas yr~'] 


kms" 


kms~ 


[pc] 


[°] 


[°] 


HD 24431 


09 IV-V 


binary 


7.2 


6.9 


-0.21±1.01 


-1.52+0.86 


-9.9+2 


12.0 


911+57 


262+38 


56 


HD 48099 


07 V 


binary 


6.3 


6.4 


+0.81±0.65 


+2.35+0.53 


+31+2 


37.7 


1829+77 


71+15 


55 


HD 57061 


09 II 


quintuple 


6.0 


6.1 


-1.82±0.44 


+3.74+0.59 


+40.4+2 


55.8 


1914+79 


116.0+6.5 


46 


HD 92206 


06? 


binary 


8.3 


8.2 


-10.4±4.2 


+6.8+4.2 


-10+5 


40.5 


4049+109 


147+19 


1 


HD 135240 


07 III-V 


triple 


5.0 


5.1 


-2.02±0.51 


-4.08+0.55 


+9.2+2 


32.5 


1131+63 


243.7+6.5 


16 


HD 149757 


09 V 


single 


2.6 


2.6 


+ 13.07+0.86 


+25.44+0.71 


-15+10 


26.8 


163+27 


62.8+1.6 


34 


HD 158186 


09.5 V 


binary 


7.0 


7.0 


+ 1.36+1.11 


-1.26+0.48 


-9 


13.9 


1094+62 


317+26 


41 


HD 17505 


06.5 V 


binary? 


7.8 


7.4 


-1.38+1.16 


-0.69+1.08 


-17+5 


25.7 


1416+69 


207+41 


41 



Table 1. Surveys containing the candidates of the Ha-sample; 
bow-shock detections as results from the overlay are printed in 
bold letters. 



Star 


Survey 


Bow Shock 


1 1 IJ 1 jj / 


V 1 M 


ambiguous 


im 1 AC 


VTCC 

V laa 


complicated background 


I 1 YJ lyoZX) 


VTCC 
V I oo 


weak 


rll) Z44jl 


v laa 


detection 


I 1 YJ JUDIH- 




not in survey 


I 1 IJ jH-U / o 


VTC C 1 


not in survey 


DD -TJ7 LJZ.O 


VTCC 
V 1 


non— detection 


i I r~v OTflOn 

rlD j /UzU 


VTC C iCU A C C A 


ambiguous 


i ir~\ a i i f. i 

I 1 YJ H 1 1 D 1 


/VTC C 1 
( V 1 o 3 J 


not m survey 


HD 41997 


VTSS 


ambiguous 


HD 47839 


VTSS 


non-detection 


HD 48099 


VTSS/SHASSA 


complicated background 


HD 52533 


SHASSA 


non-detection 


HD 54662 


VTSS/SHASSA 


ambiguous 


HD 57061 


SHASSA 


detection 


HD 64315 


SHASSA 


ambiguous 


HD66811 


SHASSA 


non-detection 


HD 92206 


SHASSA 


small 


HD 101131 


SHASSA 


non-detection 


HD 112244 


SHASSA 


non-detection 


HD 130298 


SHASSA 


ambiguous 


HD 135240 


SHASSA 


detection 


HD 329905 


SHASSA 


non-detection 


HD 149757 


SHASSA 


detection 


HD 156212 


SHASSA 


ambiguous 


HD 158186 


SHASSA 


complicated background 


HD 164492 


SHASSA 


non-detection 


HD 169582 


SHASSA 


ambiguous 


HD 175514 


VTSS 


weak 


HD 186980 


(VTSS) 


not in survey 


HD 188001 


VTSS 


non-detection 


HD 227018 


(VTSS) 


not in survey 


HD 195592 


(VTSS) 


not in survey 


HD 199579 


VTSS 


incomplete scan 


HD 203064 


VTSS 


non-detection 


HD 210839 


(VTSS) 


not in survey 


HD 214680 


VTSS 


non-detection 



Due to the high parallax errors determined by Hipparcos for 
the candidates the distances were determined using their spec- 
tral parallax with absolute magnitudes derived from Landolt- 
Bornstein ( 1982 > according to the spectral classification de- 
scribed in appendix|A] As for the absorption along the line of 
sight, photometric data for the B and V filters were taken, and 
the normal extinction law (Mathis. fl990l l applied. The expected 
absorption within the Ha-line was estimated according to the 
interstellar extinction given by Mathis ( 1990l. 
In appendixlAlthe magnitudes and multiplicity of the different 
sources are described. All data and results are given in Table|2] 

3.3. Motion 

The proper motion of the eight bow-shock candidates was 
taken from the Hipparcos-catalogue, as well as their errors de- 
rived from the given error-ellipse. The radial velocity infor- 
mation was taken from CDS (Evans D. S.[T373| Wilson R. E. 
I1953 i). Only in the case of HD 158186 was this value updated 
with respect to the ones given in VB. The astrometric data led 
to the determination of the inclination and the position angle 
of the proper motion concerning the central stars of the bow- 
shock candidates. The astrometric results are given in Tabled- 
If the nebula were created by a bow shock, these parameters 
should be the same as for the bow shocks. The inclination can 
be directly compared. As for possible image rotation of the Ha 
image compared to the global coordinate system, neighbouring 
stars were measured and the position angle corrected. 

4. Analysis 

Bow-shock nebulae are structures which appear limb bright- 
ened, due to their shell like geometry. Assuming that the gas 
within the layer is optically thin, one can compute qualitatively 
the characteristics of a brightness profile by determining the 
length of the line of sight within the layer. When plotting this 
profile against the radial distance of the line of sight, one ob- 
tains radial brightness plots as shown in Fig.|4] Case A demon- 
strates the situation of limb darkening for a sphere, while B 
shows the radial brightness plot of a spherical shell. The ra- 
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Fig. 4. Theoretical radial brightness plot for a sphere (A) light 
grey, a spherical shell (B) dark grey and a bow-shock layer 
(C) medium grey. They are derived from the line of sight s(x) 
through the regions, with the star at the position x = 0.1 and II 
are the exterior and interior of the different regions. 

HP 57061 (positional angle: 200 o -210°) 
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Fig. 5. The measured radial brightness plot of HD 57061 with 
the radial distance from the central O-star. 



dial brightness plot of a bow shock is that of case C. The in- 
ner boundary of the layer was described by a parabola, as van 
Buren et al. (1990 1 suggested. The outer boundary was de- 
scribed by a confocal parabola ensuring a constant thickness 
of the layer. The similar appearance of radial brightness plot B 
and C will be discussed later. 

For all eight candidates, detected using IRAS 60 /mi excess 
maps (see Sect. 13. U , radial brightness plots were derived by 
transforming the pixel coordinates of the images into a polar 
coordinate system centred upon the central O-star. The surface 
Ha brightness of all pixels within a 10° wide wedge were then 
plotted against their radial distance in arcmin. A resulting radial 
brightness plot for HD 57061 is shown in Fig. [5] demonstrat- 
ing the typical limb brightening of a bow shock superimposed 
upon a nearly constant background emission. 
The measurement of the symmetry axis of the bow-shock neb- 
ula was verified with a method, other than that proposed by VB. 
When applying their method one determines the symmetry axis 
of the structure created by the bow shock and the background 



emission together. This would alter the direction of the bow- 
shocks symmetry axis. If one could determine the location of 
the inner boundary of the bow-shock layer free of background 
contaminations, it would be possible to determine the symme- 
try axis of the bow shock by itself. This can be achieved using 
the location of the maximum of the radial brightness plot. The 
maximum is determined by the varying brightness of the bow 
shock alone. Background nebulosity will most certainly not be 
so sharply peaked. Therefore, its location exactly traces the in- 
ner boundary. 

The radial distance of the inner boundary was traced with ra- 
dial brightness plots for each bow-shock candidate. Plotting 
the distance against the position angle of the wedge used to cre- 
ate the radial brightness plot, results in a radial distance profile 
as shown on the right of Fig.[6]for the example of HD 158186. 
The radial distance profile shows a symmetric behaviour and 
its symmetry axis coincides with that of the bow-shock struc- 
ture. Taking the existence of a symmetry axis for granted the 
position of the axis is determined by a symmetric function fit- 
ted to the data points, which also coincides with the symmetry 
axis of the bow shock. A parabola: 

r = a{6 - 6 S ) 2 + b (1) 

was chosen to fit the data because of its few free parameters and 
its close approximation of the data. The parabola is fitted to the 
data plotted in polar coordinates. A parabola in cartesian coor- 
dinates, for which the parabola of van Buren et al. ( 1990 1 in Eq. 
|2]is defined, can be transformed, but would be of a more com- 
plicated structure not needed to measure the symmetry axis. 
Resulting from the great width of the wedge, only few radial 
brightness plots were created, thus only few data points were 
present in the radial distance profiles, but enough to ensure a 
symmetric distribution (see open or filled circles on the right of 
Fig-E). To double the data points and gain a more precise sym- 
metry axis position the points were mirrored with an assumed 
symmetry axis 8 S . Thereafter, the parabola (Eq.[0 was fitted to 
the points, using the least-square method and error weighting 
(as for all following fits), keeping the displacement 6 S constant. 
In steps of (El, G s was changed and an variance s 1 of the fit de- 
rived. The symmetry axis was chosen to be at a 6 S of minimal 
variance (see the left of Fig.[6]for the example of HD 158186), 
which agrees with the less precise value found without mirror- 
ing. The errors of 6 S result from a final fit allowing for 6 S to 
vary. All symmetry axes are given in Table [3] 
The comparison of all position angles determined either by VB 
or through radial distance profile or astrometry are given in 
Table|4] also noting the deviation of radial distance profile with 
respect to the astrometric results. The data is transformed back 
to a cartesian coordinate system with the x-axis as symmetry 
axis. Now, the physically motivated parabola of van Buren et 
al. ©90 1 



y 



-Ro 



(2) 



3i?o 

can be fitted to the data of = [-90°, +90°] nearest the apex. 
The parameter Ro so determined is known as the standoff dis- 
tance, and is used to calculate ISM densities. Eq.[2]is only cor- 
rect for an inclination of t = and thus, Rq represents an upper 
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Fig. 6. Left: The variance s 2 of different symmetry axis S referring to an arbitrary axis for HD 158186. Right: The radial distance 
profile of HD 158186 with the position angle 9 as used in the left figure. Symmetry axis is chosen as the minimum of s 2 . Open 
circles indicate the mirrored data. 



Table 3. Derived parameters of the analysed bow shocks. 8 S and 
i are the position and inclination of the symmetry axis using the 
Ha images. The standoff distances SOD are given for the two 
cases of no inclination and inclination, d is the thickness of the 
bow-shock layer perpendicular to its surface. Uncertain values 
are given without errors. 



star 




SOD(i = 0°) 


SOD(i) 


( 


d 




n 


['] 


['] 


[°] 


[10 15 cm] 


HD 57061 


202.8±2.0 


13.56+0.16 


6+11 


75+40 


16.9+4.4 


HD 92206 


79.9±4.1 


3.11+0.16 


3.9+1.8 


0+10 


9.8+2.7 1 


HD 158186 


238.3±1.7 


5.96+0.13 


7.4+1.4 





1.80+0.11 


HD 135240 


296.7±1.7 


10.65+0.16 


13.3+1.6 





2.90+0.42 


HD 149757 


67.1+1.7 


9.42+0.18 


0+48 


90 




HD 17505 


358.1+3.8 


5.82+0.18 


5.8+2.1 





0.59+0.07 


HD 24431 


312+3.1 


7.10+0.32 


5.6+13 


71+49 


6.6+1.0 


HD 48099 


238.6±5.3 


4.25+0.25 





90+0.8 


2.25+0.62 



Layer thickness of HD 92206 is given in 10 9 cm 



of the layer can be approximated by a spherical shell of thick- 
ness d. In the case of photoionisation equilibrium, 



= r 2 a HaA (T)n]d 



r'a HaA (T)EM ± (d,r) (4) 
is valid. Q is the Lyman continuum flux of the O-star as given 



for a specific spectral type by Panagia (1973) and ana a — 
5.83 x 1(L 14 cm 3 s _1 is the recombination coefficient for the Ha 
line at T = 10000K derived from Osterbrock ( TH^ . Here, 
the emission measure EM ± is only valid when looking directly 
through the bow-shock layer. The maximum EM we measure 

is given as EM(d, r) — 2EM ± (d, r) Jl^ + 1- For the parabola 
of Eq.|5]the radial distance can be calculated as: 



r(x) = 



(5) 



limit. The rotated parabola (z(x, y) = ^ - # ) with a differ- g ™ [ 6 pc cm ] 
ent inclination angle t is more complex: 



y I 3^o 7 

z(y) = cos £ — Rq I tan ( + 1 1 cos l 

3Ro \ 4 



(3) 



Both parabolas were fitted to the data of all eight bow-shock 
candidates. The standoff distance with or without inclination 
and the inclination t are given in Table |3. 
A further important parameter characterising a bow shock is 
the thickness d of the bow-shock layer. To determine d, one 
plots the variation of the brightness of the inner boundary, de- 
rived from the radial brightness plots, against the position angle 
relative to the symmetry axis resulting in a brightness profile, 
shown for HD 57601 in Fig.0 The brightness is given as an 
emission measure which can be approximated as EM = n 2 s in 
case of a homogeneous density n s within the bow-shock layer 
and the length of the line of sight s through the layer. Near the 
apex of the parabola and at an inclination of t = 0°, a segment 




Fig. 7. The brightness profile of HD 57061 resulting in a thick- 
ness d = (16.9 + 4.4) x 10 15 cm. 

and the polar angle = arctan - + 90° is given by: 



x{6) = -atanf 



-a 2 tan 2 6 + 3a 2 
2 



(6) 
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Table 4. The symmetry axis determined with the IRAS 60 pm 
excess maps (6* p ), Ha-images (# s ), and derived from the astro- 
metric data (6 a ) are given and compared using the deviation 
AO = a - S . All angles are given in deg. 



star 


e P 


6s 








A9 


HD 24431 


160 


226.0i 


3.1 


2624 


38 


36±41 


HD 48099 


340 


301.4±5.3 


71 + 15 


230±20 


HD 57061 




348.2±2.0 


116.04 


6.5 


232.2+8.5 


HD 92206 


166 


112.1±4.1 


1474 


19 


35±23 


HD 135240 


145 


250.3i 


1.7 


243.74 


:6.5 


6.3±8.2 


HD 149757 


290 


115.94 


.1.7 


62.8±1.6 


52.9±3.3 


HD 158186 




306.74 


1.7 


3174 


26 


10+28 


HD 17505 


154 


183.94 


3.8 


2074 


41 


23+45 



This together leads to a brightness profile of: 



EM(d, 0) 



+ 1 



2a HaA nr 2 (8) 



(7) 



that has to be fitted to the points of the brightness profile, after 
one has subtracted a typical value of the emission measure of 
the background measured at the edge of the images. In addition, 
an absorption correction has to be applied as well as a transfor- 
mation of angular to linear distances using the distances of the 
O-stars, both calculated in Sect. 5. The function of Eq.0was 
fitted and the results are given in Table 0. In the case of HD 
149757, the determination of d was impossible due to the sat- 
uration of the O-star contaminating the edge of the bow shock 
(see Fig.|8}. 

5. Results 

It can be seen in Eq.|5J that the standoff distance Rq is the main 
parameter determining the structure of a bow shock. The dis- 
tance of the bow-shock layer at the apex of the parabola, which 
is the standoff distance, is determined in the frame of reference 
of the central star by the balance of the ram pressure of the stel- 
lar wind and that of the moving ISM. The exact formula of Ro 



(8) 



is given by Wilkin ( 1996 



Ro = 



\47r« lSOD , y t 2 



To further derive the density of the surrounding ISM tisoDfi, 
the mass loss rate M w and the asymptotic velocity V w were 
taken from the literature (see footnotes in Table If non 
were present, the relations of van Buren ( 1983) and van Buren 
(1985 1 based on the luminosity and effective temperature as 
given by Panagia ( 1973 1 for the appropriate spectral class were 
used to calculate the missing parameters. 
Due to the large errors of the standoff distance when taking 
the inclination into account, only uninclined standoff distances 
were used to derive the ISM density. Therefore, all derived ISM 
densities summarised in Tableware lower limits. 
The temperature of the surrounding ISM is derived due to the 
fact that a bow shock can only be created when the central 
star has supersonic velocity \M = ^ > 1} with respect to the 



Table 5. The ISM parameters «sod,o and T max as derived from 
the standoff distance (Ro), compared to nu,o as derived from 
bow-shock layer density (n s ), using mass loss rates and stellar 
wind velocities taken from the literature as cited in the foot- 
notes. 



star 



[pc] 



nsoD,o 
[cm" 3 ] 



[10 3 K] [10 3 cnr 3 ] [cm- 3 ] 



HD24431 ad 


1.89: 


=0.05 


3.3= 


=0.2 


6.37 


3.914 


:0.34 


2.93±0.70 


HD 48099 bd 


2.26= 


=0.10 


0.18= 


t0.02 


62.9 


4.84 


:2.1 


1.5+1.0 


HD57061 cd 


7.56= 


=0.13 


0.07= 


=0.002 


138 


4.284 


:0.76 


12.4±5.4 


HD 92206 aa 


3.67= 


=0.19 


0.007= 


=0.001 


7260 


4810+340 


5.2+1.8 1 


HD 135240 bd 


3.50: 


=0.08 


0.21 = 


=0.01 


46.8 


4.564 


-.0.95 


1.43+0.50 


HD 149757 bd 


0.45±0.02 


1.5= 


=0.1 


31.8 








HD 158186 aa 


1.90: 


=0.05 


2.0= 


=0.1 


8.55 


7.34 


:3.2 


2.96±1.5 


HD 17505 bd 


1.20= 


=0.04 


21 = 


tl 


29.2 


23.64 


-.1.9 


2.47±0.50 



1 ISM density of HD 92206 is given in 10~ 3 cnT 3 

a van Buren ( 1983 l and van Buren ( 1985 l 

b Howarth & Prinia fE%9l 

c Lamers & Leitherer <1993t 

d Howarth et al. (19971 

surrounding ISM. The local sound velocity within a neutral 
medium is 



C,= 



ykT 



(9) 



whereby k is the Boltzmann-coefficient and the adiabatic ex- 
ponent is assumed as y — |, because the mean mass per parti- 
cle and proton mass is /j. — 0.61, when the Helium fraction is 
0.1. As any velocity higher than sound speed can create a bow 
shock, the derived temperature is only a maximal one and also 
given in Table [5] 

As one has derived the thickness d, measured the maximum 
EM of the bow-shock layer, and therefore EM ± , the density of 
the layer can be derived using the definition of the EM= n s d. 
The surface density at the apex (cf. Wilkin, 1996) 



3 / V, 

cr = n s d = -Ran LT ,o 1 1 + — 



(10) 



then leads to an alternative calculation of the ISM density Hup 
(Table|5}. 



6. Discussion 

The results are discussed in three parts. First, the individual 
characteristics of the different candidates are mentioned and 
the problems encountered while analysing them. Second, the 
quality of the developed methods to detect bow shocks and de- 
termine their parameters are discussed. And thirdly, the sample 
and its characteristics as a whole are described and elucidated. 



6.1. Candidates 

As can be seen in Fig. [5] and Table 0] the symmetry axis of 
the bow-shock candidate around HD 48099 deviates greatly 
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Fig. 8. Median filtered Ha images of the first four bow-shock candidates. The short arrow indicates the proper motion direction 
of the star and the long arrow the position of the symmetry axis. The position of the star is given by a star symbol. The grey-scales 
used are given in rayleigh (R) or decirayleigh (dR). The inset shows the IRAS 60 /im excess images of the bow-shock candidates 
used to determine the overlays, their size scaled to the Ho- images and using inverted grey-scales. 



from the proper motion direction. Due to the neighbouring ISM 
cloud which dominates the bow-shock nebula in the northeast- 
ern side, the axis is located along the main axis of the neigh- 
bouring ISM cloud. 

In case of HD 149757 the saturation of the central star causes 
'bleeding' through the bow-shock nebula (see Fig. [8} and leads 
to an erroneous direction of the symmetry axis (see Table|4}. As 
before the deviation is in the expected direction, opposite to the 
defective region. 

For HD 57061 no explanation can be found for the deviation 
of the symmetry axis. The nebula shows typical properties of a 
bow-shock nebula like its limb brightening (Fig.|5}, as well as 
the double hump feature at the apex, seen in Fig. [8] predicted by 
Mac Low et al. ( 1991 1 for bow shocks inclined like HD 57061. 



This could be a hint to a possible misdirection of the astromet- 
rically determined proper-motion direction of HD 57061. 
The bow shock around HD 92206 could not be resolved com- 
pletely. Therefore the derived ISM parameters deviate from the 
expected values for the Warm Ionised Medium (WIM). 
Comparing the symmetry axis 9 S from the Ha images with 
those derived from the astrometric data 6 a in Table |4]it can be 
said, that they coincide within their errors. Only the above men- 
tioned cases show a significant misalignment. The same can be 
said, when comparing the ISM parameters with expected val- 
ues for the WIM. All in all the case of the eight candidates 
being bow shocks can be significantly strengthened. 
As for the spatial velocity noted in Table |3 the following can 
be said: the large errors of the velocities of HD 149757 and 
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HD 17505 
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Fig. 9. Median filtered Ha images of the last four bow-shock 
defined as in Fig. [8] 

HD17505 are sufficient to classify them as runaways. As for 
HD 149757, it is a long known runaway star (cf. Hoogerwerfer 
et al. 1200 1> . The determined spacial velocities of HD 24431 
and HD 158186 are much lower than the runaway limit, thus 
the runaway character is uncertain. However the nebulae de- 
tected are most certainly bow-shock nebulae. The remaining 
four candidates have typical velocities for a runaway OB-star. 
The radial brightness plot derived using the Ha images only 
shows the expected profile of a bow shock in case of a few 
position angles. All profiles only show the brightening of the 
nebula toward its limb, due to the strong background emission. 
These profiles are sufficient to determine the position of the in- 
ner boundary of the bow-shock nebula, but are useless as a sole 
indicator for a bow shock. 

The best indicator is the correlation of proper-motion and 
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candidates with their IRAS 60/mi excess images, symbols are 



symmetry-axis direction. Deriving radial distance profiles and 
their symmetry-axis have been shown to be robust against 
bright background emission, through the use of radial bright- 
ness plot, as mentioned above. The deviations from the ex- 
pected profile can easily be analysed and put into context with 
the surrounding ISM, for example when encountering an ISM 
gradient. 

Comparing ISM parameters derived from the bow-shock struc- 
ture with expected WIM parameters is also a promising method 
of clarifying the bow-shock character of the nebula. Only the 
inclination-free fit of the bow shock leads to a stable fit. The 
determination of the standoff distance with respect to the incli- 
nation leads to faulty results, as all candidates seem to either 
have inclination of nearly 0° or 90°. Such a high fraction of 
highly-inclined and non-inclined bow shocks is quite improb- 
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Fig. 10. The stellar distance histogram of the complete Ha sam- 
ple. The bow-shock detections are plotted in black. 
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Fig. 11. The stellar space velocity histogram of the complete 
Ha sample, as in Fig. 1101 



able. Especially as the inclination within the used sample is 
not only uniformly distributed, but biased by a preselection of 
visible bow shocks, such will have low inclinations improving 
their visibility. 

Deriving ISM parameters using the determined thickness of the 
bow-shock layer and comparing them with the WIM, is also a 
good method of clarifying a nebula as a bow shock. The thick- 
ness is more difficult to determine, being not directly measur- 
able. It is derived by fitting an appropriate brightness profile, 
which results from a crude model of the brightness distribution 
along the bow-shock layer. The density of the layer has been 
taken as constant, which is only correct in the vicinity of the 



apex, see Wilkin (1996 1. Due to the changing angle between 
bow-shock layer and ISM ram pressure, the force bounding 
the layer will also change, resulting in a change of thickness. 
As the brightness profile fit was only done near the apex of the 
parabola, these effects are not as dominant as the errors made 
when flux correcting. The correction of the background emis- 
sion is done by approximating a constant emission, leading to 
great uncertainties of the derived thickness. However the ISM 
parameters derived are within their errors comparable to the 



(1 = 1 80) 




A bowshock 



© 0B association 

HI shells, superbubbles 



Fig. 12. The distribution of the bow shocks in context to super- 
bubbles and OB-associations, as published in Huthoff & Kaper 
(2003. 



ISM parameters derived from the standoff distance. 
The remaining problem and limit of all methods is their ability 
to discern bow shocks from stellar bubbles. As demonstrated 
in Fig.|4]the radial brightness plot shows a clear limb brighten- 
ing in both cases. Only the difference in the overall geometry 
(axi-symmetric in the former or spherically symmetric in the 
latter) enables us to find the bow shock. But, when taking the 
surrounding ISM having a density gradient, the structure of a 



stellar bubble described by Castor et al. ( 1975) is altered. The 
resulting stellar bubble becomes ellipsoidal and the part at the 
high-density end will be brighter or the only visible part. Such 
a deformed stellar bubble can hardly be discerned from a bow 
shock, except for its symmetry axis lying in the direction of the 
density gradient. The most noticeable difference is the veloc- 
ity of the matter in the layer for both cases. In the case of a 
stellar bubble, the matter moves radially away from the central 
star, whereas the matter in the bow-shock layer moves along 
the layer and approximately tangential to the star. Hence, only 
velocity information of the matter within the layer, gained ei- 
ther by spectroscopy or velocity charts within H i can solve this 
last ambiguity (see e.g. Brown & Bomans 2003 1. 

6.2. Sample 

The distance histogram containing all eight detections is com- 
pared to that containing 22 non-detections, without taking into 
account multiplicity and absorption. The resulting histogram, 
shown in Fig.^| reveals that most detections lie within a dis- 
tance of 2 kpc. This cannot be a result of surface brightness, as 
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it would stay constant with respect to varying distances. The 
bow shock, or rather the standoff distance, will vary in size 
when the distance changes. For example HD 92206 is the most 
distant detection and hence is hardly resolved. Similar results 
are gained when analysing the velocity histogram in Fig. II llfor 
the sample. Only bow shocks around slow central stars can be 
resolved. One can say this sample is complete within 2 kpc and 
up to velocities of 60 km s _I . 

As described in appendix|A] the sample contains only one sin- 
gle system. As for the DES and BSS, only 50% of the cre- 
ated runaways should be multiple systems. Whether more bow 
shocks are created by multiple systems as for single star sys- 
tems, or if the fractions of multiples determined from both sce- 
narios is wrong, cannot be stated on the basis of only eight bow 
shocks. 

The location of the eight bow shocks within the galactic plane 
in Fig.^]shows that they are not located inside any superbub- 
ble. This is to be expected, as the sonic velocity (Eq. |9] with 
T ~ 10 6 K) therein would be to high to create a bow shock. 
Only HD 17505 seems to lie in the superbubble inside GS 137— 
27-17, which is only an effect of projection. HD 17505 having 
a b — -0.9° lies well above the upper surface of the bubble at 
b = 0°. For the three bow shocks with 270° < I < 0° no data 
for superbubbles exists. 

The typical lifetime of an OB-star with its runaway velocity 
does not permit it to put a large distance between itself and 
the OB-association as a probable origin. As shown in Fig. 1121 
most bow shocks are still near OB-associations, except HD 
92206, HD 57061 and HD 17505. The first two, are the fastest 
stars within the sample, enabling them to move a greater dis- 
tance, and HD 92206 is also so far away, that OB-associations 
can only be determined with great difficulty at such a distance. 
HD 17505 is in the vicinity of a star formation region noted in 
Carpenter et al. (2000 1 not shown in Fig. 1121 

7. Conclusion 

The search for bow shocks using the Ha-sample of SHASSA 
and VTSS yielded eight detections (HD 17505, HD 24430, HD 
48099, HD 57061, HD 92206, HD 135240, HD 149757, and 
HD 158186) from a total of 30 candidates already observed 
(seven candidates are missing in VTSS), derived from the sam- 
ple VB used for their search within the IRAS allsky survey. 
The best indicator to detect bow shocks within Ha images was 
the correlation between the direction of the proper-motion and 
the symmetry axis, determined using radial distance profiles, 
which are not sensitive to bright backgrounds. The other meth- 
ods can be used for further verification, when in doubt (as for 
HD 57061). 

The detected bow shocks could be successfully used to deter- 
mine ISM parameters. This was done, either using the stand- 
off distance or the brightness profile. As the sample is only 
complete up to a distance of 2 kpc and no bow shocks can be 
found inside a superbubble, the derived values of the density 
(~ 1 cirT 3 ) and the maximal temperature (~ 10 4 K) fit well to 
the picture of the WIM (e.g. Shull ll987> . Both features jus- 
tify the more or less constant ISM density and temperature. 
Regions of other ISM composition and therefore other param- 



eters, e.g. the regions inside the Small and Large Magellanic 
Cloud, are too far away and cannot be analysed using the 
medium resolution of the Ha-allsky surveys. All in all, bow 
shocks around OB-runaway stars are ideal probes of the ISM, 
as are neutron stars (cf. Chatterje & Cordes 2002). They also 
demonstrate the qualitative picture of the neighbouring ISM is 
apparently correct. 

First steps to clarify the last ambiguity of bubble or bow shock 
using spectroscopy or velocity charts (see Sect. 16. 1> were taken 
for the eight candidates: Using data from the IUE-archive, ab- 
sorption profiles of N v, Si iv, and C iv were detected in the 
case of HD 48099, HD 57060, and HD 135240 showing an ex- 
cess compared to surrounding O-stars and low temperatures 
of ~15 000K. As stellar winds are too thin to contribute large 
amounts of these ions and the temperatures are too low to ra- 
diatively excite them, this confirms again these objects as shock 
fronts of bow shocks. More important, these lines were used to 
measure the velocities, which fit to the proposed values of the 
bow-shock layer. Additionally, Hi velocity data of HD 17505 
from the Canadian Galactic Plane Survey was used to create ve- 
locity charts, which showed a velocity distribution as expected 
in the case of a bow-shock layer. However no data could be 
found in the literature concerning the other candidates. Hence 
further spectroscopic analysis of these objects is still needed 
to make use of the full potential of the current sample as ISM 
probes. 
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Appendix A: Photometric Data 

HD 135240: Penny et al. J200T1 have determined HD 135240 
as a triple system and specified the first component as a 07 
III-V, the second as a 09.5 V and the last one as a B0. 5 V 
star. They also measured the UV flux ratios for the different 
components as ^ = 0.239 and = 0.179, so that the 
magnitudes of the single components could be determined 
from the total magnitudes given in the SIMBAD-database. 
HD 57061: For photometric data, the SIMBAD-database was 
used and the results of Stickland et al. J1996i of the star being 
a quintuple system. They could resolve one component as a 
09 II star. The other consists of two double systems resulting 
in a total spectral class of B0. 5 V. The single O-star is the 
dominant star with a ten-times higher flux than the second 
double binary. 

HD 17505: Fabricius & Makarov (2000) note this star as being 
a double or triple system and could measure the magnitudes 
of the brightest component using the Tycho-filtersystem. The 
spectral class of 09.5 V (Garmany et al. 1982) can only be 
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determined for the complete system. 

HD 158186: This object was detected as being a variable 
of the Algol-type by Adelmann et al. (2000 1, therefore HD 
158186 has to be at least a double system. As the system is not 
resolved, the star was treated as a single star with photometric 
data from the SIMBAD-database and the total spectral class 
of 09.5 V given by Buscombe ( fT998l . 



HD 24431: Fabricius & Makarov (2000) also analysed this 
star and classified it at least as a double system and were able 
to measure the magnitude of two components. As in the case 
of HD 17505 they used the Tycho-filtersystem. Reed {T998i 
states the spectral class as a 09 IV- V. 

HD 92206: The SIMBAD-database notes this star as a double 
system, but can only give the total magnitude of the system. 
Therefore it is treated as a single star. Reed ( 1998 1 could only 
measure the spectral class of 06, and was not able to determine 
the luminosity class. 

HD 149757: Here the SIMBAD-database notes the star as 
a single system and its magnitude. The spectral class was 



measured by Garmany ( 1982 1 and is 09 V. Being the nearest 
star of the eight candidates, the distance determined by 
Hipparcos could be determined. Both distances are within 
reasonable agreement with respect to their errors, verifying the 
spectral parallax results. For consistency, we use the distance 
determined by spectral parallax. 

HD 48099: This star is given as a binary system by Stickland 
et al. (1996 1. They could measure the UV-flux ratios for 
the components as = 1.8. Garmany i 19821 1 could only 
determine the total spectral class to 07 V. 
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